Atoms in double-5-kicked periodic potentials: chaos with long-range correlations 
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We report an experimental and theoretical study of the dynamics of cold atoms subjected to 
closely-spaced pairs of pulses in an optical lattice. The experiments show the interplay between fully 
coherent quantum dynamics and a novel momentum-diffusion regime: for all previously-studied 5- 
kicked systems, chaotic classical dynamics shows diffusion with short-time (2 or 3-kick) correlations; 
here, chaotic diffusion combines with new types of long-ranged 'global' correlations, between all 
kick-pairs, which control transport through trapping regions in phase-space. Analytical formulae 
are presented and, with quantum localization, are used to analyse the experiments. 

PACS numbers: 32.80.Pj, 05.45.Mt, 05.60.-k 



The '<5-kicked particle' (<5-KP) is one of the most stud- 
ied experimental and theoretical paradigms of classical 
Hamiltonian chaos. A particle is kicked periodically by a 
sinusoidal potential V(x,t) = —KcosxJ2^(t ~ n )- F° r 
sufficiently large K, the motion is fully chaotic and the 
energy grows diffusively. Of particular interest in recent 
years has been the theoretical pj, |2| and experimental 
[3j investigation of the suppression of the classical dif- 
fusive process in its quantum counterpart, the quantum 
8- kicked particle (<5-QKP). This phenomenon is gener- 
ally termed 'Dynamical Localization' (DL). For K >^ 1 
the diffusive growth of the classical energy is no longer 
bounded by phase-space barriers (tori) so persists indef- 
initely : for an ensemble of classical particles, the energy 
E =< ^- >= \Dt ~ ^-t for all t. The corresponding 
quantum energy grows only up to a timescale t* ~ D/h 2 
0| and saturates at a value < P 2 > t _ +00 ~ Dt*. The 
Texas experiments showed that cesium atoms in period- 
ically pulsed waves of light were an ideal test-bed for 
quantum chaos. A broad range of interesting physical 
regimes were subsequently investigated: controlled deco- 
herence Q, quantum accelerator modes and dereal- 
ization induced by non-periodic kicking Q. 

The chaotic diffusion is however not entirely uncorre- 
lated ja| and there are corrections which have now been 
experimentally observed due to 2-kick and 3-kick cor- 
relations. For example, a 2-kick correction appears be- 
cause the ensemble averaged value for the correlation be- 
tween the impulse at the n — th kick and that experienced 
2-kicks later Co =< V' (x n )V' (x„.-i-2) > is generally non- 
zero, see 8, 121. In hfi hj it was further shown that, 
if the pulses are unequally spaced, the 2-kick corrections 
yield a local (in momentum) correction to the diffusion. 
In that case, 2E(Po,t) = D(Po,t) t : for unequal kick 
spacings, both the diffusion rate and hence the energy, 
depend non-trivially on time and the relative initial mo- 
mentum, Pq, between the atoms and the standing wave 
of light. 

In this work we report the first experimental and the- 
oretical study of the 25- kicked particle (2<5-KP): a cloud 
of cesium atoms is exposed to a periodic sequence of 
closely spaced pairs of kicks. At the outset, one might ex- 
pect that the diffusive behaviour here could be analysed 



within the framework used in |gL 1 1 0l| - of diffusion with cor- 
relations between short sequences of 2 or 3 kicks, whether 
local or otherwise. However this approach fails to explain 
the experimental results. Further investigation showed 
that chaotic diffusion in the 25-KP was rather different 
from that seen in previously studied kicked systems. For 
the 2<5-KP one finds new corrections, which appear in 
families correlating all kicks democratically. These cor- 
rections are individually very weak, but become numer- 
ous with time and accumulate to eventually dominate 
the diffusive process. Moreover, these 'global' correla- 
tions can be associated with specific physical phenomena, 
namely the escape from and through 'trapping' regions 
in phase-space. We have identified diffusive regimes asso- 
ciated with three types of correlation: one is an ordinary 
1-kick correlation, C± , the other two are new and are fam- 
ilies with 'global' terms. The experimental behaviour de- 
pends on which correlations dominate at the point when 
the quantum dynamics suppresses the diffusion. 
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FIG. 1: Surface of Section diagram for an ensemble of par- 
ticles, with K = 7 and e = 0.1, prepared with initial mo- 
mentum Po = at t = 0. The SOS illustrates the trapping 
of trajectories in phase space regions for which the momenta 
pe ~ ±(2n + 1)tt. 

The experimental apparatus is essentially like that de- 
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scribed in |ll| and consists of a cloud of cesium atoms 
collected in a standard 6-beam MOT and cooled in an op- 
tical molasses to a temperature of 6/zK. The sinusoidal 
potential V(x,t) is formed by two counter-propagating 
laser beams incident on the cloud with parallel polarisa- 
tions. These are derived from a Titanium Sapphire laser, 
have an intensity of 4 x 10 3 P a t (I sa t =1.12mWcm~ 2 , the 
saturation intensity) in each beam and are detuned 2000r 
(r = 2vr x 5.22MHz, the natural linewidth) below the D2 
transition on cesium. The potential is switched on us- 
ing acousto-optic modulators (AOMs) to create pulses as 
short as t p = 300ns, and each beam is controlled by a 
separate AOM so that a frequency difference A/ may be 
imposed upon the two beams and the potential moves 
with constant velocity in the laboratory frame. In this 
way we may explore the momentum dependence of the 
diffusion constant as in the rest frame of the potential 
the atomic momentum distribution has a non-zero mean 
value P oc A/. 

We now have two periods: r, which represents a (long) 
timescale between the pairs of kicks, and e, which repre- 
sents a much smaller time interval between kicks in each 
pair. In the experiment, r = 9.47/iS, while five different 
separations of the closely-spaced pair, in the range 0.44/zs 
to 1.48/iS were used: e = 0.047, 0.063, 0.094, 0.125 and 
0.156. For these parameters and with the intensity and 
detuning as above, we have an effective value of K = 1 
while the kick-strength, K = 3.3 (±10% due mainly to 
the uncertainty in measuring the intensity in the laser 
beams). Up to 100 kicks were applied before the cloud 
of atoms was allowed to evolve freely in the dark for 
15ms. A pair of near-resonant imaging beams were then 
switched on and the fluorescence imaged on a CCD cam- 
era. From the spatial distribution of the fluorescence the 
momentum distribution was extracted and (P 2 ) calcu- 
lated. 

The corresponding classical behaviour of the 2<5-KP 
would be given by evolving a 2-kick map: 

Pi = Pi-i + Ksinxi-i 

Xi = Xi-l +Pi€ 

Pi+i = Pi + Ksinxi 

X i+ 1 = Xi +p i+ lT 

It is instructive to begin by considering the classical evo- 
lution of an ensemble of particles all initially at time t = 
0, with momentum p = Po for which Pqc = (2n+ l)ir and 
n = 0, 1, ... These particles experience a kick K sinx^-i 
followed by another at ~ K sm(xi^i + tt) which in effect 
cancels the first. Conversely in the case Pq6 = 2mr a 
series of near-identical kicks produces initially rapid en- 
ergy growth. Other Po produce intermediate behaviour. 
This behaviour follows from the unsurprising fact that 
the presence of the short time interval, e, results in a 
non-zero kick-to-kick correlation C\. In contrast, for the 
Standard Map and all other previously studied atomic 
<5-kicked systems, C\ = 0. 



We show, in FigQ] Surface of Section plots obtained 
from an ensemble intially with Po — 0, but randomly 
distributed in position x, for K = 7, e = 0.1. It is clear 
that, though the phase-space here is fully chaotic, tra- 
jectories 'stick' at the values of pe ~ ±(2n + 1)tt. By 
calculating G\ explicitly, after N pairs of kicks, one can 
begin a more precise analysis: 

Ci(Ar,P ) = A- 2 cosP e[Jo(A-e) - J 2 {Ke)\ 

N 

Y^iMKe))^- 2 (0.1) 

3=1 

Physical time is t = N(t + e). We can re-scale variables 
so r = 1 and t ~ N. The ensemble- averaged energy E of 
the atom cloud at time i, corrected by Eg 10. II would be 
given by 2P =< (P t -P ) 2 >~ K 2 T/2 + d(t, P ) (where 
T = 2N = 2t). It is easily shown that for short times C\ 
grows linearly in time, while for longer times it saturates 
to a constant value, after a time t\ ~ r^jp ■ O n this 
time-scale, the kick-to-kick contributions decay to zero. 
Since C\ is a single correlation term, analogous to those 
studied in [10| we have here been able to simply extend 
the usual analysis of 0,0]. For small Ke, J 2 (Ke) ~ 0, 
so for short times (t << ti) we can write for the energy, 
2E =< (P t - P ) 2 >~ K 2 T /2[1 + cosP e]. The validity 
of this formula for £ << ti is seen in the experimental 
data in FigE^a). 

FigElshow the energy absorbed by the cesium atoms, 
as a function of Po. Three separate values of e were 
considered. For each e, the energy was measured, after 
100 kicks (but here t* ~ 40) for many values of Po. In 
every case, at time t = 0, the atom cloud had energy 
< (Pt=o — Po) 2 >— 30, indicated by the horizontal dot- 
dashed line. 

The experimental data of Figl^a) is the most straight- 
forward to understand: here, t\ ~ 450, so t* << t\ 
so the energy absorption was arrested in the regime 
t << t\. An ensemble of classical particles, initially at 
momenta pje ~ ±(2n+ l)ir, absorbs no energy while for 
p)£ ~ ±2rt7r, energy absorption is maximal. 

FigEIc) , on the other hand, was the most surprising. 
It corresponds to a regime t* > t\. It shows a clear 
reversal of the behaviour seen in Fig|2Ia): atoms initially 
prepared at or near the momentum-trapping regions end 
up with more energy than those prepared in the enhanced 
momentum diffusion regions. In fact, atoms which are 
prepared furthest from the momentum trapping regions, 
absorb the least energy. 

This counter-intuitive observation motivated a more 
careful look at the mathematical model of the dif- 
fusive process of the 2<5-KP. This exposed a collec- 
tion of terms involving averaging products of the form 
sin Xi sin x^ , where fi < i but is otherwise arbitrary; 
while i represents any of the second kicks in the 
pair. For each [i we obtain a correction of the form 
4if 2 cosP e Jl{Ke)Y, j {Ja{Ke)) 2j - z . These terms are 
negligibly small (O(Ke) 2 relative to C±). However since 
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FIG. 2: Experimental results for 25-KP realisation with ce- 
sium atoms. Each data point (star) shows the energy ab- 
sorbed (after 100 kicks, K = 3.3, h = 1) by a cloud of atoms 
with average momentum p = Po (relative to the optical lat- 
tice) at initial time, t = 0. With increasing e, we see the 
minima (maxima) in the energy flip into maxima (minima) 
as the 'global' correlation family Cgi gradually overtakes the 
nearest-neighbour correlation Ci. The dashed lines represent 
a classical simulation using 100,000 particles all with mo- 
menta = Po at t — 0, and K within the range 3.3 ± 10%. 
(a) t* << ti ~ l/(Ke) 2 . Regime dominated by the one-kick 
correlation C\. Atoms prepared near the trapping regions 
(Poe ~ (2n + l)7r) remain trapped. Results follow closely the 
formula < P 2 >~ K 2 T/2(1 + C osP e). (b) t* ~ l/(Ke) 2 . 
Regime for which Ci competes with the 'global' correlation 
family Cgi- This begins to expose the inverted peaks of the 
Poisson correlation terms Cp, which determine the trapping 
very close to the resonant condition (P e = (2n + l)7r). (c) 
t* > l/(Ke) 2 . Re gime dominated by Cgi , but sharp inverted 
peaks due to Cp are still visible, (d) Competition between 
Ci, Cgi and leading order Poisson terms. At short times, 
the global correlations are negligible, but at later times, the 
accumulation of (individually weak) global correlations over- 
take the corresponding Ci and Cp m terms: when the global 
terms become dominant, all the atoms have escaped from the 
trapping regions. 



we sum over all /it < i their numbers accumulate with 
time and the net contribution of this 'global' correlation 
family is: 

- C G1 (N,P ) = 4K 2 cosP e J 2 (Ke) 
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It is easily shown that, though negligible at short times, 



FIG. 3: The figure plots the energy absorbed by an ensem- 
ble of 100,000 classical particles, all with p = Po at t = 0, 
as a function of initial momentum (dashed line). These nu- 
merics are superposed with analytical formulae obtained for 
the diffusion (solid line). K = 7 and e = 0.05 in every fig- 
ure, but the number of kicks, 2N, varies. Note the behaviour 
in the trapping region Po ~ n/e ~ 60: the energy absorbed 
is a minimum at short times but turns into a maximum at 
long times, (a) Results at short time (10 kicks), dominated 
by the one-kick correlation C\. (b) Results at a time for 
which Ci competes with the 'global' de-phasing family Cgi- 
This exposes the 'Poisson-sum-formula' corrections: a fam- 
ily of corrections of the form Cp m oc cos mP e, (m > 2) 
combine to give a series of inverted peaks in the trapping 
regions (c) All three corrections Cp, C\ and Cgi compete (d) 
At t > 500 the longest lived correction, Cgi dominates. The 
Poisson corrections have changed sign: the inverted peaks of 
(c) (indicating trapping), have 'flipped over' to give a series of 
positive peaks which give the maxima a 'pointed' appearance. 
The counter-intuitive result, that particles initially prepared 
in the momentum-trapping regions will be the ones which will 
acquire the most energy, is evident. 



this term grows quadratically at small t and eventually 
overtakes C\. It has opposite sign to C\ ; we interpret it as 
a term which reflects the gradual de-phasing of the reso- 
nant effects of C\ -such as the kick-cancellation/trapping 
regions with Pqe ~ (2n + l)7r). At long times the be- 
haviour is dominated by Cgi and hence we see that the 
energy absorption is maximal for particles prepared near 
P e ~ ±(2n+ 1)tt. 

Fig. 2(b) corresponds to a particularly interesting 
regime, where C\ <~ Cgi- The cosPoe correction is 
partly cancelled, exposing a series of narrow dips in 
the energy. The origin of these dips is in a series of 
terms Cp m oc cos mP e. When summed these pro- 
duce behaviour reminiscent of the Poisson sum formula 
E m (- 1 ) mcos ^oe = E n *(^oe-(2n+l)^). The ampli- 
tudes of the C p m terms vary with time and only a finite 
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number of harmonics (m < 10 typically) contributes at 
any given time. Hence we get a series of broadened peaks. 
Nevertheless, for this reason, we term these corrections 
the Poisson term, Cp = ^2 C'p m and where: 

m— 1 

C Prn (N, P ) ocK 2 cos mP Q eF m (t) [] jf(nKe). (0.3) 

n=l 

F m (t) is a time function which grows as ~ t m . Though 
these terms are small (0(Ke) 2m ~ 2 ) they will contribute 
when t m (Ke) 2m ~ 2 ~ 1. Each of these terms also 
has a partner 'global' family of opposite sign, Cgf™ 
(0(Ke) 2m ), which increase as ~ t m+1 . Fig. 2 (d) shows 
the behaviour of the m = 2 Poisson terms in compari- 
son to C\ and Cgi. As m increases terms become less 
significant and for each m the global correlation always 
dominates Cp m at long times. There are additional, even 
higher-order, cosnPoe terms (n > 1) of similar form to 
Cp m and CcPm above and a group of terms involving 
products of the form sin 2 x% which may contribute, par- 
ticularly to Cp2- 

In Fig|3|we compare numerical (classical) calculations 
at K = 7, e= 0.05, with the behaviour predicted by the 
correlations. In Fig[3Ja) we look at short times and see 
that we can accurately match the energy absorption by 
including only the three lowest order diffusive corrections, 
ie those which increase linearly or quadratically in time 
(Ci, Cgi, Cp%). At later times (t < 50 or so), we obtain 
good agreement by including all the above terms up to 
order Jf° and m = 4: the inverted peaks correspond- 
ing to the trapping regions are quite well reproduced, as 
seen in Fig|3Jb). At extremely long times, we can obtain 
good results simply from the leading global family Cgi- 
We also clearly see, in Fig^d), the inversion of the be- 
haviour seen at earlier times: atoms prepared within the 
trapping regions give positive peaks, since CcPm domi- 
nates. The — th order term is < P 2 >~ K 2 T/2 at very 
short times, then gradually slows down to an asymptotic 
value ~ K 2 T/16 (obtained numerically) at long times. 
The analytical curves in Fig. 3 are shifted vertically by a 
constant amount. 

In conclusion, we have presented and analysed an ex- 
perimental realisation of a 26 — KP and shown that it 
corresponds to a type of diffusion quite different from 
the Standard Map. This is a rich and complex system 
and many questions remain open. On the classical side, 



certain aspects of the interactions and lifetimes of the 
different types of correlations are not yet well under- 
stood. Although we have a generic handle on this sys- 
tem with the diffusion correlations, further insight might 
be gained by detailed knowledge of classical trajectories. 
For lower values of Ke, a thin band of islands and eventu- 
ally unbroken tori appears first around pe ~ ±(2n+ 1)tt. 
At this point, clearly, there will be no escape through 
the trapping regions and beyond the very shortest times 
(ii << l/(Ke) 2 ) the diffusive correlation approach would 
fail everywhere. We expect regions permeated by broken 
phase-space barriers to persist at the parameters such as 
eg. Fig^and to account for the trapping regions. 

The quantum behaviour is also not fully understood. 
Though clearly some key features can be modelled qual- 
itatively in Fig. 2 by adjusting the number of kicks and 
K (by up to 10% in the classical numerics of Fig. 2), the 
quantum dynamics will involve tunnelling and dynamical 
localization (DL) effects which can introduce substantial 
differences with the classical behaviour. Further investi- 
gation of the DL in this system is necessary: its Floquct 
states are localised, like in the standard S — KP. But 
in the latter, localization lengths, L, are quite uniform 
with exponential localization and L ~ while in the 
28 — KP, values of L can vary by about 3 orders of mag- 
nitudes for the parameters considered here |T3 |. 

Although the t « ti experimental regime (seen in 
Figs. 2(a) and 3(a)) is in some sense the least surprising 
in terms of the diffusive process, it is worth noting its 
potential for atomic manipulation. Interest in the 2<5-KP 
experiment was initially stimulated by its potential ap- 
plications in manipulating atoms in devices like an atom 
'chip' |l3|. In |l(J it was proposed that a local diffusion 
rate D(P Q ), could be exploited for 'filtering' cold atoms 
according to their velocity. For selected Pq the atoms 
could pass the device unperturbed, while other momenta 
would absorb a substantial amount of energy and would 
be dispersed. The 2<5-QKP showed a far stronger experi- 
mental signature than the system in fldi Hll | which relied 
on a two-kick, C%, correlation (note that C2 7^ Cpi). 
A much stronger velocity-selective effect, due to the Ci, 
correlation is seen for the 25-QKP. The inverted peaks of 
the C p correlations could also be used to select a narrow 
band of velocities with Pq ~ ir/e. 

The authors thank Thibaut Jonckheere for helpful dis- 
cussions. This work was supported by the EPSRC. 
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